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ABSTRACT 

We study the dynamical origin of the structures observed in the scattered- 
light images of the resolved debris disk around HD 141569A. The disk has two 
conspicuous spiral rings and two large-scale spiral arms. We explore the roles of 
radiation pressure from the central star, gas drag from the gas disk, and the tidal 
forces from two nearby stars in creating and maintaining these structures. The 
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disk’s color, scattering function, and infrared emission suggest that submicron¬ 
sized grains dominate the dust population observed in scattered light. CO obser¬ 
vations indicate the presence of up to 60 of gas. The dust grains are subject 
to the competing effects of expulsive radiation pressure (/? > 1, where f3 is the 
ratio of the radiation and gravitational forces) and retentive gas drag. We use 
a simple one-dimensional axisymmetric model to show that the presence of the 
gas helps conhne the dust and that a broad ring of dust is produced if a central 
hole exists in the disk. This model also suggests that the disk is in a transient, 
excited dynamical state, as the observed dust creation rate applied over the age 
of the star is inconsistent with submillimeter mass measurements. We model in 
two dimensions the effects of a fly-by encounter between the disk and a binary 
star in a prograde, parabolic, coplanar orbit. We track the spatial distribution 
of the disk’s gas, planetesimals, and dust. We conclude that the surface density 
distribution reflects the planetesimal distribution for a wide range of parameters. 
Our most viable model features a disk of initial radius 400 AU, a gas mass of 
SOM®, and /? = 4 and suggests that the system is being observed within 4000 yr 
of the fly-by periastron. The model reproduces some features of HD 141569A’s 
disk, such as a broad single ring and large spiral arms, but it does not reproduce 
the observed multiple spiral rings or disk asymmetries nor the observed clearing 
in the inner disk. For the latter, we consider the effect of a 5 Mj^p planet in an 
eccentric orbit on the planetesimal distribution of HD 141569A. 

Subject headings: hydrodynamics — planetary systems: formation — planetary 
systems: protoplanetary disks — stars: individual (HD 141569) — circumstellar 
matter 


1. Introduction 

Debris disks around main-sequence stars are dusty, optically thin, and gas-poor. Radi¬ 
ation pressure (RP) and Poynting-Robertson (PR) drag eliminate dust grains on timescales 
shorter than the stellar age, so the observed dust must be continuously replenished by col¬ 
lisions among, or evaporation of, planetesimals (Hackman & Paresce 1993). The Infrared 
Astronomical Satellite (IRAS) and the Infrared Space Observatory (ISO) revealed over 100 
stars with far-infrared excesses indicating the presence of debris disks. However, spatially 
resolved images of these disks are relatively rare: only about a dozen debris disks have been 
resolved since the early 1980s (see Zuckerman 2001 and references therein). 

The resolved disks are not featureless. They frequently display warps, spiral structures. 
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and other azimuthal and radial asymmetries. Different mechanisms have been suggested to 
explain these features. They may be caused by interactions between the dust and the gas 
in the disk (Takeuchi & Artymowicz 2001, hereafter TAOl), the formation of small planets 
(Kenyon & Bromley 2004), the dynamical forces of embedded planets (Ozernoy et al. 2000), 
or stellar-mass companions (Larwood & Kalas 2001; Augereau & Papaloizou 2004; Quillen 
et al. 2005). Whatever their cause(s), the features observed in the dust disks provide insight 
into the characteristics of the unseen planetesimal population and illuminate the dynamical 
processes of young planetary systems. 

HD 141569A (AO V; age 5±3 Myr, Weinberger et al. 2000) has a resolved circumstellar 
debris disk. Its Hipparcos distance is 99 pc from the Sun and has co-moving M2 V and 
M4 V companions separated by 1.4” and located at distances of 7'.'55 and 8793, respectively 
(Weinberger et al. 2000; Augereau et al. 1999). The disk was hrst resolved in scattered light 
by Weinberger et al. (1999) and Augereau et al. (1999) using the Hubble Space Telescope’s 
(HST’s) Near-Infrared Camera Multi-Object Spectrometer (NICMOS). Mouillet et al. (2001) 
imaged the disk with HST’s Space Telescope Imaging Spectrograph (STIS) and noticed 
strong brightness asymmetries. These observations have recently been complemented by 
images from HST’s Advanced Camera for Surveys (ACS) (Clampin et al. 2003, hereafter 
C03). Together, these scattered-light images reveal very complex structure, including an 
inner clearing within 175 AU of the star, a bright spiral “ring” with a sharp inner edge from 
175 to 215 AU, a faint zone from 215 to 300 AU, and a broad spiral “ring” from 300 to 
400 AU (C03). These distances are measured along the projected disk’s southern semimajor 
axis. (For convenience, we will hereafter refer to the two tightly-wound spiral rings simply 
as “rings.”) C03 also observed two low-intensity, large-scale spiral arms in the outermost 
part of the disk, that they attributed to tidal interaction with the M dwarf companions. 

The influence of the M dwarf companions on the disk’s morphology was first discussed 
by Weinberger et al. (2000). They argued that if all three stars were coplanar with the 
disk, the system would not be stable. Given its young age, however, the system may be 
bound. With this assumption, Weinberger et al. (2000) concluded that resonant interactions 
between the companions and the disk do not account for the disk structure. Assuming that 
the companions revolve around the primary star in a highly eccentric orbit, Augereau & 
Papaloizou (2004) constructed dynamical models of the optical disk comprising dust grains 
that respond to the stellar gravitational held. They assumed that the grains are large enough 
that RP, PR drag, and gas drag are unimportant. They also ignored the presence of the large- 
scale spiral arms reported by C03. Their model reproduces the general appearance of the 
brightness asymmetries observed in the disk. Quillen et al. (2005) constructed hydrodynamic 
models of HD 141569A’s gas disk without dust, and concluded also that the companions he 
in an eccentric orbit. Their models reproduced the large-scale spiral arms (even after several 
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periastron passages) and some other disk asymmetries. 

TAOl modeled the interaction between disks of gas and intermediately-sized (> 8 /rm) 
dust grains, and Lecavelier Des Etangs, Vidal-Madjar, & Ferlet 1998 considered disks with 
very small particles. It is not clear that the conditions explored by any of these models 
is the correct one for the scattered light images of the disk around HD 141569A. In §2, 
we argue that submicron-sized grains account for most of the scattering opacity at optical 
wavelengths. Such small grains are subject to RP which quickly expels them from the disk. 
The gas detected in the disk (Zuckerman et ah 1995) tempers somewhat the short “blow 
out” timescales, but its mass is not enough to dominate the dynamics of the dust (§3). In 
this paper, we explore the effects of strong RP, gas drag, and gravity on small dust grains, 
using HD 141569A’s disk as a benchmark. We simultaneously track the behavior of three 
different populations - the unseen planetesimals controlled by gravity, the gas controlled by 
gravity and volumetric fluid forces, and the dust controlled by gravity, RP, and gas drag. 
The insights obtained from this analysis should be applicable to other systems where RP is 
important. 

In §2, we review the observational constraints that every model of the system should 
satisfy. We consider evidence that indicates the dominance of submicron-sized grains and we 
set limits on the amount of gas in the disk. In §3, we explore the interaction between gas and 
dust by means of a one-dimensional model. We show that the gas slows the outward motion 
of the dust and can produce broad ring-like structures. In §4, we present two-dimensional 
dynamical simulations of the interactions between planetesimals, dust grains, gas, and the 
three stars. Assuming that the two M dwarf companions are co-moving but unbound, we 
find that a recent parabolic fly-by causes some of the observed disk structure, including the 
large-scale arms. However, our models produce a more disorganized disk than is observed. 
We argue that the central hole in the disk can be produced by a planet with a mass a few 
times that of Jupiter (or, alternatively, some number of smaller planets) in a highly eccentric 
orbit about the primary star. 


2. Observational Constraints 
2.1. Grain Size 

Figure 1 shows the optical depth prohle of the disk, taken from Figure 5 of C03. This 
prohle represents the median values of concentric annuli centered on the primary star. Drawn 
from a composite of F435W (ACS B band) and F606W (ACS broad-R band) images, it is 
roughly indicative of the optical characteristics of the disk at 0.5 /rm. Power-law hts to 
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the inner and outer edges of the prohle (using as errors the standard deviation at every 
radius) give for r < 200 AU and for 320 < r < 500 AU. Beyond 500 AU, the 

average azimuthal brightness seems to be dominated by the light from HD 141569BC. The 
optical depth prohle indicates that the amount of mass in ~ 0.5 /xm-sized grains is ~ O.OIM 0 
(assuming a constant dust opacity of 2 10^ cm^/gr and albedo of ~ 0.5. See Wood et al. 
2001 ). 

The disk is occulted by ACS’s coronagraphic mask within ~ 150 AU. Fisher et al. (2000) 
and Marsh et al. (2002) detected mid-infrared thermal emission within this radius and Marsh 
et al. (2002) concluded that the optical depth at 1.1 /xm decreases by a factor of ~ 4 within 
this region. Li & Lunine (2003, hereafter LL03) modeled the thermal spectrum of the disk 
by assuming a density prohle similar to the one shown in Figure 1. Their model suggests 
that the reduced emission within ~150 AU is due to reduced dust density and not a change 
in the scattering properties of the grains. 

C03 reported that the disk is redder than the star, with color excesses of A{B—V) = 0.21 
and A(U — /) = 0.25. They also reported no color variation as a function of distance from 
the star (which supports our use of a constant opacity in the mass calculation above). C03 
inferred that the disk’s colors are consistent with astronomical-silicate grains having a size 
distribution of and a minimum radius of s ~ 0.4 /xm. Augereau & Papaloizou (2004) 
used the same color information to derive lower size limits between ~ 0.1 /rm and 3.1 /rm. 
The existence of such small grains is also implied by the small scattering asymmetry factor 
(0.15 < g < 0.25^) of the Henyey-Greenstein function derived by C03. The multicomponent 
models by LL03 suggest that the minimum grain size is between 0.1 /rm and 10 /rm: the 
smaller limit produces too little emission in the IRAS 60 /rm band and the larger limit 
produces too much emission in the IRAS 60 /rm and 100 fim bands. Taken together, these 
arguments imply that there is a population of grains whose radii may extend down to 0.1 fim. 

The combined influences of the reduced scattering efficiency {Qsca) of small grains and 
the small number of large grains (assuming a size distribution going as produce sharply 
peaked scattering opacity, which indicates that the grain sizes responsible for the optical 
images of the disk are of the order of the wavelength of observation. Figure 2 shows the 
scattering opacity at 0.5 fim as a function of grain radius for astronomical silicate grains (Laor 
& Draine 1993; Draine & Lee 1984) larger than 0.01 fim. The scattering opacity is expressed 
as the scattering cross section per unit mass (oc ascaS~^, where o-gca = QscaTTS^), weighted 
by mass and number of grains at each radius s (Miyake & Nakagawa 1993). The function 
is centered at 0.2 fim with a characteristic width of As ~ 0.3 fim. The sharp decrease in 


^Because of a typo, the value of g for the HD 141569A disk is quoted as g = 0.25 — 0.35 in COS 
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the scattering opacity of grains larger than 0.2 /xm is determined by the behavior of the first 
scattering peak in Qsca- In the limit of large grains, the scattering opacity decreases more 
slowly for intervals of constant ds /s). For this dust model, 1 yum grains scatter ~ 8 

times less than 0.2 /im grains, and ~ 3 times more than 10 /xm grains. The location and 
width of the peak in the scattering opacity depends on the surface properties of the grains 
and the exact mixture of silicates, ices, and empty space. For example, increasing porosity 
moves the peak to smaller sizes and less-reflecting grains broaden the peak and move it to 
larger sizes (Bohren & Huffman 1983). 

If - as suggested by the disk colors, the phase of the scattering function and the spectral 
energy distribution - submicron-sized grains exist in the disk, their dynamics are strongly 
influenced by RP. According to TAOl, the ratio {(3) of the forces from RP and gravity is ~ 4 
for a 1 fim dust grain, and /? ~ 20 for s ~ 0.05 /rm. For the very porous grains assumed by 
LL03, /3 ~ 12 for s ~ 1 fim. Formally, and in the absence of gas, grains with f3 > 0.5 (“/5- 
meteoroids”) are expelled from the disk. Grains with large (3 will acquire a radial velocity 
on the order of their Keplerian velocity in Tesc ~ (-P/27r)//5, where P is the orbital period at 
a given position (TAOl). At the edge of a 500 AU disk around a 2.3 Mq star, Tesc ~ 300 yr 
for = 4. Such short expulsion timescales are somewhat mitigated by the effects of another 
force: the gas drag. 


2.2. Gas Content 

Zuckerman et ah (1995) measured the gas content of HD 141569A’s disk using the 
detected J = 2 —1 transition for and the upper limit on the same line for ^^CO. 

They assumed that the star’s distance was 200 pc, its disk had a radius of 130 AU, and that 
H2/CO ~ 10000. They calculated gas masses between 20 and 460 for optically thin and 
optically thick limits, respectively. Assuming that the gas and dust disks are coincident, we 
recompute the gas content using the observed extent of the optical disk and the measured 
Hipparcos distance to HD 141569A of 99 pc. 

Zuckerman et ah (1995) have generously provided us with their unpublished digitized 
CO spectra, which are reproduced in Figure 3. The spectra show a double peak characteristic 
of sharp disks, although the depths of the central line-emission in both measurements are a 
little over la per resolution element. Assuming an disk inclination of 55° (C03), the velocity 
difference between the peaks is 3.8± 1.0 km s“^, which implies emission at 380^jgQ AU. (The 
upper limit is dictated by the size of the optical disk.) Assuming H2/CO ~ 4000 (Lacy et al. 
1994), the optical depth of the central line is r = 1.5 x 10® (S/g^r cm“^) (5.5A'/T)®/^, where 
S is the surface density and T is the excitation temperature of the gas (Beckwith & Sargent 
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1993). Nondetection of the line yields an upper limit of T = 30 mK, so the opacity of 
the line is uncertain. If the line is optically thick, then we obtain a gas mass of 

< 60 M 0 . If the line is optically thin and hlls a disk of radius 380 AU, we obtain a gas mass 
of 0.02 M 0 . In each case, the gas is assumed to be in local thermal equilibrium (LTE). 

3. The Dynamics of High /? Particles 

We now explore the interaction between small dust grains, gas, and the radiation and 
gravitational fields of HD 141569A, while ignoring the presence of the binary companions. 
Most of the results of this axisymmetric treatment should be applicable to other debris and 
gas disks in which RP is important. The regime that we explore here is that of modest 
amounts of gas (up to 50 M 0 ) and not very large values of /3 (/3 ~ 4 — 10). Our calculations 
are complementary to those of Lecavelier Des Etangs et al. (1998), who consider (3 > 100 
and TAOl, who consider /5 < 1. Notice that Lecavelier Des Etangs et al. (1996) has shown 
that if the observed dust particles all have (3 < 0.5, a single ring of planetesimals in a gasless 
disk will produce a dust disk with surface density given by r“^, as observed here. 

We solve the equations of motion for a dust grain, assuming that it is subject to RP, 
PR and gas drag. The gas is assumed to be in a circular orbit and we assume the same dust 
model as TAOl. The force of the gas on the dust is given by 

Fg =Av, (1) 

where Pg is the gas density, s is a dust grain radius, is 4/3 times the thermal velocity of 
the gas, and Av is the relative velocity of the dust with respect to the gas. To determine 
the trajectory of a dust grain, we assume that the surface density of the gas (as well as 
that of planetesimals) decreases radially as This function is characteristic of the dust 

surface density of optically thick protoplanetary disks (Osterloh & Beckwith 1995). In other 
words,we assume that the gas and planetesimal density prohles of the optically thin disks are 
the same as the dust density prohle of the optically thick disks. We further assume (unless 
otherwise stated) that the gas is in LTE, with temperature 

Tg = 278 K (2) 

where L* = 22 . 4 L 0 is the stellar luminosity (Backman & Paresce 1993). This implies that 
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Without gas, the radial velocity of a dust grain is given by 
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where Tq is the radius at which the dust grain is created and Vo is the Keplerian velocity at 
To- In general, there is not an analytic expression for the radial velocity of a dust grain in the 
presence of gas. For the range of parameters we are considering here, the dust grain speeds 
are supersonic and almost radial at the outer edge of the disk (Figure 4). The dimensionless 
stopping time (the stopping time in units of the inverse local Keplerian frequency) at the 
edge of the disk is given by (TAOl): 


ApdSVK 

1 S ^ - 

SpgVVr 




-)( 


gas 


Vr/Vk' 


( 5 ) 


which for the parameters considered here is < 1. Notice that, unlike in the case considered 
by Lecavelier Des Etangs et ah (1998), the radial velocity shown here does not have a 
deceleration region, because in our case the gas surface density decreases with radius: the 
radiation pressure always dominates the dynamics of the dust. 

The effect of gas on the grains depends on where the grains were created. Figure 5 
shows that for a hxed value of (3 and without gas, particles that start closer to the star soon 
overtake particles farther out. This changes with the presence of gas, because the drag force 
decreases outward, and for > 5 OM 0 masses of gas, the particles created closer in will not 
overtake those created farther out. Figure 5 also shows that the gas mass, not the grain size, 
determines the escape timescale: there is not much difference between (3 = A and /? = 10 for 
50 Mq of gas, because both the gas drag and the RP are proportional to (3. 

If one assumes a rate of dust generation, the results of these calculations can be used to 
predict the steady-state dust surface density in an isolated disk. In such a state, grains are 
continuously generated by collisions and then blown away by RP. The dust surface density 
(assuming conservation in the number of particles) is given by 


S(r) oc 



dNjro) 
Vr{r^y r)r’ 


( 6 ) 
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where r* is the inner radius of dust generation, Vo is the position at which the grains are 
created, dN{r-o) is the number of grains created at r,,, and Vrivo, r) is the radial velocity at r 
of grains created at Tq. The dust creation rate per unit volume is taken to be ^rei x o'cozz x Np^, 
where Np is the number of planetesimals per unit volume, acoii is the collision cross section 
and Yrei is the relative velocity between the planetesimals, which is proportional to the 
dimensionless planetesimal disk thickness {Hp/r) times the Keplerian velocity (Thebault 
et ah 2003). This implies 

dN{ro) = Vrei CTcoii Np^ Hp 27rrodro oc ffk Sp 27rrodro (7) 

where Sp is the surface density of planetesimals at the creation point, which we assume to 
be proportional to r~^'^. 

In practice, this simple prescription (the “quadratic” dust generation prescription) may 
be affected by uncertainties in the planetesimal size distribution and in the relative velocities 
among the the planetesimals. To explore the sensitivity of the results to the exact generation 
mechanism we also consider a “linear” dust generation prescription in which dN{ro) oc 
'EpTodVo. 

For these calculations we assume that the planetesimals and the gas are well mixed and 
therefore Hp = H. The resulting steady-state surface densities are shown in Figures 6 and 
7. In the hgures, the unit length is 1000 AU. Figure 6 shows that even modest amounts of 
gas have a signihcant effect on the surface-density prohles. The amounts of gas considered 
(10 and 50 M^) are not enough to conhne the dust grains created at distances larger than 
0.01, but gas drag sharpens the prohles, as it slows down the dust particles. Notice that if 
dust creation starts very close to the central star, the resultant surface-density prohles are 
ehectively featureless. Figures 6 and 7 also show that a broad ring of dust is produced by 
truncating the dust creation within a certain radius. Just outside this limit, the number 
and surface density of planetesimals, and hence dust grains, are large. Far from the creation 
limit, the surface density of the dust decreases because of decreasing surface density of 
planetesimals and the geometric dilution of “blown out” dust grains created at smaller radii. 
Thus, if a mechanism exists for preventing the creation of dust at small radii (for example, 
a planet that clears out the parent planetesimals), the result is, very naturally, a broad ring 
of dust at larger radii. 

The rate at which dust is lost from the disk depends on the rate of planetesimal erosion. 
In principle, one can use use the former to estimate the latter. LL03 calculated the amount 
of mass lost from HD 141569A’s disk, assuming no gas, a grain-size distribution oc (3=1 
for grains of all sizes, and a constant rate of mass loss throughout the age of the system. 
They also assumed a surface-density prohle for the dust that is slightly different than the one 
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shown in Figure 1. They concluded that 39 M0 of solid material have been lost, in particles 
with 1 < s < 10 /xm, from the disk due to RP and PR drag. 

To examine the effect of the gas on the rate of mass loss, we consider a disk of radius 
500 AU about HD 141569A and no binary companions. Figure 8 shows the time required 
for small dust grains to move from a given radius to the edge of the disk for different values 
of gas mass and temperature. The presence of just 10 M0 of gas triples the escape time 
from 200 AU, and a colder disk conhnes the dust more than a hotter one. (Equation 1 shows 
that, neglecting ur, a colder ~ denser - gas disk, produces more drag). Assuming a grain-size 
distribution oc with 50 M0 of gas in LTE, varying escape times, the TAOl dust model 
and a constant rate of mass-loss, we determine that ~ 700 of solid material have been 
lost over the ~ 5 Myr age of the disk, in particles with sizes 0.1 < s < 8 /rm. Without gas, 
~ 1200 M0 of solid material would have been lost. The difference with LL03 is due to the 
fact that in our calculation different size particles have different blow-out timescales. 

However, neither LL03’s estimate nor ours can be correct. Because small grains (< 
10 /im) are quickly expelled from the disk, they must be the dominant component of the lost 
mass. If the grain-size distribution is oc for grains with radii < 1 mm, then grains up 
to 10 /im compose only ten percent of the disk mass. The bulk of the dust mass comprises 
large grains that are eliminated in much longer timescales. If the collisional processes re¬ 
sponsible for the replenishment of small grains also produce large grains, then submillimeter 
observations should detect ~ 10^ of dust, which greatly exceeds the amount of 2 M0 
measured by Sylvester, Dunkin, & Barlow 2001. 

One possible solution to this discrepancy is that the dust is not generated with a power- 
law size distribution, but with a distribution weighted toward small particles. The collisional- 
evolution models of Thebault, Augereau, & Beust 2003 suggest that, although large varia¬ 
tions in the dust distribution are possible, most of the mass produced by the collisions is in 
large particles. These models assume that the collisional timescale is smaller than any other 
timescale in the disk. 

Another possible solution is that the current rate of dust creation is larger than it has 
been in the past. This possibility is consistent with the idea that the disk has been recently 
stirred by a close encounter with a companion or unbound star. In section 4 we show that 
a parabolic encounter would occur over timescales on the order of 10^ yrs. Over this time 
the disk would lose (according to the model above) ~ 1 of solids in small particles, or 
~ 10 M0 in particles up to 1 mm in size, close to the measured value. Observationally, the 
presence of ~ 0.01 M0 of solids in 0.5 /rm-sized particles (Section 2) suggests that there is 
currently on the order of 0.1 in particles up to 8 /xm in size, with the exact number 
depending on the assumed value of the opacity and dust size distribution. 
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Furthermore, these estimates support the idea of a single exciting event, like an en¬ 
counter with an unbound companion, as opposed to repeated exciting events. In the models 
by Augereau & Papaloizou (2004), which assume a bound companion, repeated encounters 
are required to explain the brightness asymmetries in the debris disk. 


4. A Dynamical Simulation 

The presence of large-scale spiral arms suggests that the dynamical influence of the 
binary companions, HD 141569B and C, is important. If they are coplanar with the disk, 
their center of mass is ~ 1250 AU away from the primary star and their separation is 
~ 275 AU. Augereau & Papaloizou (2004) assumed that the companions are bound to 
HD 141569A in a very eccentric orbit. Their dynamical models show that close encounters 
between the companions and the disk produced well-developed spiral rings within the disk 
after only a few periastron passages. However, the models do not maintain large-scale spiral 
arms after a few orbits. Thus, within their model, the interaction geometry needed to 
create the spiral rings in the disk is not consistent with the presence of the large-scale spiral 
arms. The hydrodynamic model by Quillen et al. (2005) does produce spiral arms after 
repeated close encounters, because the viscosity and gas pressure help generate a spiral at 
every periastron passage. However, the amount of gas actually present in the disk (§2.2) 
discourages the notion that the dynamics of the dust conforms to that of the gas. Figure 
5 shows that, for a reasonable range of gas mass, /?-meteoroid grains responsible for the 
scattered-light disk are impeded by gas drag, but they are not bound by the gas. So the 
assumption by Quillen et al. (2005) that gas dynamics control the observed structures in the 
dust disk is not appropriate. 


4.1. Assumptions and Methodology 

We assume that the encounter between the disk and the binary companions is a parabolic 
fly-by, a situation known to produce spiral arms in dusty disks (Larwood & Kalas 2001). 
Such an encounter is consistent with the proper motion of the system. Within an error 
box of O'.'l, the position of the companions relative to HD 141569A has not changed over 
the 60 yr baseline of observations noted by Weinberger et al. (2000). The largest expected 
relative motion is obtained when the system is currently at periastron and the encounter 
occurs in the plane of the sky. In this scenario, the relative motion of the companions over 
60 yr would be 07012, which is well inside the error box of the proper motion measurements. 



Dust grains are created by collisions of planetesimals. We assume that the planetesimals 
themselves experience gravity, but not RP or gas drag. We also ignore self-gravity and the 
dynamical effects of the collisions. Our simulation is performed in two phases. First, we 
obtain the distribution of planetesimals and gas as a function of time. Then, for each 
planetesimal configuration, we generate a dust distribution and follow it as a function of 
time. We repeat this sequence for all planetesimal and gas conhgurations. At any given 
time, the dust in the disk consists of dust created at that time and dust remaining from 
earlier times. 

The model is two-dimensional. We assign hxed masses of 2.3, 0.5, and 0.25 Mq to 
HD 141569A, B, and C, respectively. We assume that the encounter is prograde, as a 
retrograde encounter fails to produce large-scale spiral arms. The phase of the encounter is 
hxed so that, at periastron, the three stars are aligned, with HD 141569B between A and 
C. (The results are mostly insensitive to the choice of phase.) Most of our simulations start 
with 10^ planetesimals. The inner dust destruction limit is set at r > 0.1, where the distance 
r is measured relative to the separation of HD 141569A and HD 141569BC’s barycenter at 
periastron. We assume that the initial surface density of planetesimals declines as The 

outer radius of the planetesimal disk is a free parameter. Figure 9 shows the evolution of 
the planetesimals in a disk of unit radius. We dehne T = 0 as the time of periastron. The 
time unit is such that the period at unit distance is ) = 7.2 time units, where 

mi = 2.3 Mq and m 2 = 0.75Mq (the later is the sum of the companion masses). Notice that 
the particular fly-by conhguration shown in Figure 9 shows that some planetesimals can be 
captured by the companions and therefore that some of the IR excess measured by IRAS 
may be associated with them. 

The total mass of gas is a free parameter. We consider gas disks with r > 0.01, an 
initial surface-density prohle oc r~^'^ and masses of 0, 10, 50, and 100 M 0 , before truncation 
by the encounter. We model the gas evolution using the Smoothed Particle Hydrodynamics 
(SPH) formalism (Monaghan 1992), which allows one to track the motion of gas pseudo¬ 
particles using Lagrangian equations. Viscosity is parametrized according to Equation 4.2 
of Monaghan (1992). SPH permits the exploration of parameter space without a large 
investment of computer time, but with limited resolution at small scales. 

We explore different limits for the gas sound speed, which is parametrized by the disk 
opening angle, H/r = Cs/tVLki where H is the scale height of the gas, Cg is the speed of sound, 
and VLk is the Keplerian angular frequency at radius r. We consider gas with H/r = 0.05, 
0.1, 0.2, and in LTE (Equation 3). Brandeker et ah (2004) has shown that the empirical 
opening angle of the (I Pictoris gas disk is H/r ~ 0.28, which is close to the expected LTE 
value at 1000 AU. Figure 10 shows the evolution of gas particles for the case of H/r = 0.1. 
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To simulate the dust grains, we consider the linear (particles created with r > 0.1) 
and quadratic (particles created with r > 0.15) dust-creation prescriptions described in §3. 
Particles are created every 0.01 time units. The ratio (3, which couples the dust and the 
radiation, is also a free parameter. We explore two values: (1) /5 = 4, which corresponds to 
s ~ 1 /im (TAOl) or s ~ 4 /rm (LL03), and (2) f3 = 10, which corresponds to s ~ 0.3 /rm 
(TAOl) or s ~ 1 /xm (LL03). We assume a single value of (3 per simulation, which is 
equivalent to assuming that the scattering opacity (Figure 2) is a delta function. 

To obtain the dust distribution at a given time, we consider all the dust particles created 
at previous times. Figure 11 shows the behavior of dust particles created at various time 
intervals for /? = 4, Mgas = 50 M^, H/r = 0.1, and linear dust creation. The generations 
of particles were tracked to T = 1.5, and their respective surface densities were summed 
and azimuthally averaged at each time interval. This simulation tracks the entire history of 
dust generation, but in practice, dust blowout eliminates from the “current” surface density 
prohle all the contributions of dust created before the last half time unit (i.e., T — 0.5). 


4.2. Simulation Results 

We now present our models for the fly-by encounter with a composite disk of planetes- 
imals, gas, and dust. Representative simulations are shown in Figures 12-16. Although 
most of the models are shown at T = 1.5, the images and prohles are representative of 
0-8 ^ T < 2.0. To give our distance and time scales physical meaning, we set the distance 
between HD 141569A and the barycenter of the companions at T = 1.5 equal to that ob¬ 
served in the ACS images (C03). Thus, one distance unit corresponds to 718 AU and one 
time unit corresponds to 1800 yr. The companions are out of the frames. 

For each of the Figures 12-16 we show the predicted logarithm of the density prohle 
with linear dust generation (top row) or quadratic dust generation (bottom row). The use 
of the logarithmic stretch allow us to highlight faint features. The images are all normalized 
to the same arbitrary constant. In general, the color table of the bottom row has a lower 
limit value that the color table at the top row. This means that “white” in the bottom 
panels corresponds to somewhat less material than it does in the top panels. In this way, 
we can show faint features in the bottom rows that would not be visible were we to use 
the same color table as that from the top rows. This modihcation has little effect on the 
color the dark features, which implies that all panels within a given hgure can be compared 
to each other. As mentioned before, an artihcial inner hole has been set for the linear (at 
r = 0.1) and quadratic simulations (at r = 0.15). In each panel, the plot inset compares 
the measured density prohle with the model density prohle. The latter is normalized to the 
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value measured at 500 AU. 

Figure 12 shows the evolution of the dust structures for /5 = 4, Mgas = 50 M0, H/r = 
0.1, an initial disk size of one unit distance, and both linear and quadratic dust generation. 
Figures 13 and 14 show the effect of the amount of gas on = 4 and /? = 10 dust particles, 
respectively. The sharp, concentric rings in the gasless models (panels 1 and 4) are artifacts of 
discrete dust generation. Figure 15 illustrates the role of the gas temperature, and Figure 16 
shows the roles of disk sizes and planets. A side-by-side comparison of the observed disk and 
our favorite model is presented in Figure 17. 

All the models look similar. Despite the dynamical importance of RP, the overall ap¬ 
pearance of the dust disk resembles the underlying distribution of planetesimals. Figures 13 
and 14 show that models with gas reproduce the observed prohles better than models with¬ 
out gas, but it is not clear that 10 M0 of Hjr = 0.1 gas (panel 2 of each hgure) is better than 
50 M0 of H/r = 0.1 gas (panel 3 of each hgure). There are no signihcant qualitative differ¬ 
ences between the (3 = A and /? = 10 conditions shown in Figures 13 and 14, respectively, 
as expected from the results in §3. There is also no clearly preferred prescription for dust 
creation. The quadratic prescription does produce a more disorganized structure, as the dust 
is mostly generated in the denser spiral arms. Overdensities in the planetesimal distribution 
produce streams of matter in these models. Nevertheless, the combination of large (3 and 
quadratic dust generation with H/r = 0.1 gas (Figure 14, panels 5 and 6) reproduces very 
well the decay of the surface-density prohle. Notice that because more gas conhnes the dust 
better, the spirals in panels 3 and 6 of each hgure look darker than the others. 

Figure 15 shows that the simulations are sensitive to the gas temperature. Cold gas 
(panels 1 and 4) tends to produce sharper spiral structures than hot gas, because cold gas 
is denser and imparts greater drag on the dust. This ehect is less pronounced for quadratic 
dust generation, as the large density variations obscure the temperature-dependent ehects 
of the gas. 

All the models produce a two-armed spiral. The arms are less pronounced for quadratic 
dust generation than for the linear dust generation because the former emphasizes regions 
of large planetesimal density. The density and extent of the arm diametrically opposite the 
companions can be reduced (and made more compatible with the observations) by diminish¬ 
ing the size of the initial disk (Figure 16). The appearance of this arm can also be changed 
with a different encounter geometry. Larwood & Kalas (2001) show that, in a encounter 
between a circumstellar disk and a passing star, the spiral arm opposite the projectile is 
affected by the indirect tidal component of the potential, which diminishes as the collision 
becomes more hyperbolic. 
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In general, the models with quadratic dust generation produce more disorganized images 
than models with linear dust generation. Planetesimal condensates are strong sources of dust 
that produce “noise” in the models. However, they also produce more peaked surface-density 
prohles, which match the observed prohles better. Flatter initial density prohles (not shown 
here) produce large amounts of dust at large distances. Disks with initial radii larger than 
~1000 AU produce remnant structures that are not observed, including very long spiral arms. 
Conversely, disks with an initial radius of 400 AU yield dust distribution that reproduce well 
the observed large-scale arms and surface-density prohles (Figure 16). We surmise from 
our models that HD 141569A’s initial planetesimal and gas disks had radii < 700 AU and 
10-50 Mq of H/r = 0.1 gas (Figure 17). 


4.3. The Effect of Giant Planets 

From Figure 11 it seems clear that different dust conhgurations (i.e. large holes) can 
be created by manipulating the dust generation history of the system. The size of the hole 
can be controlled by varying (3 or the time between dust generations. Though theoretically 
possible, such contrived scenarios are physically unreasonable. 

A more reasonable agent for creating a hole in the dust distribution is a planet (or 
planets) within ~ 150 AU. If large enough, planets would clear out the planetesimals and 
the gas, and create a region without dust. Each planet would sweep up a lane along its orbit 
as wide as a few times its Hill radius, Rh- From Bryden (G. and Lin), a 5 Mjup planet with 
a semimajor axis of 100 AU and eccentricity e = 0.6 would clear a region with a half-width 
of 90 AU. Alternatively, multiple planets in less eccentric orbits would have similar effect. 

Here we present exploratory simulations that include the effect of a giant planet. Pan¬ 
els 2 and 4 of Figure 16 show models that include a 5 Mjup planet in an eccentric orbit 
(e ~ 0.6) around HD 141569A with a semi-major axis of 100 AU. The planet clears the 
central region enough to greatly reduce the dust generation close to the star and produce 
an inner hole in the dust disk. Weak signatures of the eccentricity (the elongated shape and 
slight dust overdensities at the top and the bottom of the hole) are observed. The inner 
ring in the observed density prohle is not reproduced in the observations, at least with these 
orbital parameters. Considering the uncertainties in the exact dust conhguration close to 
the coronagraphic mask, this is a satisfactory result. A 5 Mjup planet of age 5 Myr would be 
~ 10“^ times less luminous, and over 12 magnitudes fainter in 77, than HD 141569A (Burrows 
et ah 1997; Golimowski et al. 2004). Such an planet would be below the non-coronagraphic 
NICMOS detection limits (Krist et al. 1998). 
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The formation of small planets may produce sharp rings in the disk by stirring nearby 
planetesimals and enhancing dust generation. This process has been modeled in detail by 
Kenyon & Bromley (2004). Both this effect and the presence of a giant planet modeled above 
produce inner ring variability on timescales of ~ 1000 yr. 


5. Discussion 

Our favored model of the scattered-light disk is shown in Figure 17. This model assumes 
an initial disk size of 400 AU, /3 = 4, and gas with with mass 50 and an opening angle 
of Hjr = 0.1. The model reproduces the large-scale spiral arms and, despite a large value 
of /3, produces localized structures. It does not reproduce the sharp and coherent rings or 
disk asymmetries of the disk around HD 141569A, as observed by 003, which suggests that 
a different encounter geometry should be considered. The model includes a 5 Mjup planet 
in an eccentric orbit (e = 0.6) with apastron 160 AU for the purpose of clearing out a region 
of planetesimals close to the central star. The characteristics of this planet are not tightly 
constrained by the model, and other conhgurations (including multiple planets) are possible. 
The simulations clearly indicate that a mechanism capable of clearing the inside of the disk 
(like a planet or planets) is necessary to explain the observed density prohle. There may be 
additional observational evidence supporting the existence of a planet. Brittain & Rettig 
(2002) have detected within 7 AU of HD 141569A, which they suggest may come from 
the extended atmosphere from a giant protoplanet. 

Although there is favorable evidence for the existence of planet (s), we should also con¬ 
sider the effect of ice condensation on the surface density of the dust. An abrupt change in 
the surface density is expected upon crossing the “snow line,” which for HD 141569A’s disk 
is ~ 150 AU (LL03). The location of the snow line depends strongly on grain size, so it likely 
varies with observed wavelength. Current coronagraphic images lack the spectral resolution 
needed to probe these variations. While crossing the snow line cannot explain the paucity 
of dust close to the star, it likely affects the shape of the observed surface density. 

Most of the gas simulations do a similar job reproducing the optical depth images and 
prohles. It would seem that the only hard constraint that they provide is that the presence 
of gas is necessary. The quadratic dust generation method does a somewhat better job 
than the linear one in matching the density profiles. Interestingly, Figure 18 shows that 
the planetesimal distribution itself can also match the outer surface density prohle. The 
surface density of planetesimals starts as but the truncation by the encounter makes 

it steeper by the time the simulations are examined. The surface density of dust generated 
quadratically also matches the distribution of planetesimals: while the dust is generated as 
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the square of the planetesimal distribution, it is blown outward by radiation pressure. This 
process tends to flatten the simulated prohle. 

The applicability of our models to the debris disk around HD 141569A depends on the 
premise that the observed dust grains are small enough so that j3 > 1. This condition is 
borne out by the disk colors, the mean-scattering phase, and thermal models. Could larger 
particles {j3 < 0.5) be contributing to the observed scattered light image? In deriving the 
scattering-opacity function (Figure 2), we assumed that the collisional equilibrium timescale 
is smaller than any other systemic timescale, which implies that the number of grains of size 
s decreases as If small grains were more quickly removed from their points of creation, 

then an excess of larger grains would result because they would no longer suffer erosion from 
collisions with the smaller grains (Thebault et ah 2003). In the absence of gas, (3 < 0.5 grains 
are eliminated by PR drag on timescales orders of magnitude larger than the RP timescales 
(Chen & Jura 2001), so the proportion of larger grains would further increase. However, in 
the presence of modest amounts of gas, (3 < 0.5 grains migrate to stationary orbits (TAOl), 
and small particles settle down at the outer edge of the disk. The timescale over which 
this migration occurs is quite small for small grains, but larger than the purely dynamical 
timescales associated with (3 > 0.5 motions. All these processes may alter the shape of the 
dust size distribution. If, for example, the distribution is flatter than then the peak of 
the scattering-opacity function is broadened. While the optical image of the disk may still 
be dominated by /5-meteoroids, the contributions of larger, bound grains may be signihcant. 

The entangled roles of RP, gas drag, and the dynamical role of the companions greatly 
complicate any effort to model all components of the disk. To do so successfully requires 
multi-wavelength images of the disk. Observations by the Spitzer Space Telescope should 
soon help this situation, as Spitzer will observe larger grains that are less affected by RP. 
Future observations with the Atacama Large Millimeter Array (ALMA) will also contribute 
valuable information about the composition and distribution of the dust. 


6. Conclusions 

We have developed dynamical models to investigate the structures observed in the debris 
disk around HD 141569A. The models include, for the first time, the effects of radiation 
pressure and gas drag. The disk’s colors and scattering phase indicate that the scattered- 
light images are produced by submicron-sized dust grains. To understand these images, we 
must consider the behavior of grains with values of (3 (the ratio of radiation and gravitational 
forces) larger than one. These grains are expelled from the disk in shorter-than-dynamical 
timescales. We show that up to 60 M 0 of gas may be present in the disk. Gas drag slows 
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the blow-out of grains and is fundamental to understanding the structure of the disk. 

If the dust is continuously created, its steady-state surface density is featureless and has 
a radial prohle dependent upon the amount of gas and the prescription for dust creation. To 
produce a broad ring of dust, its creation must be prevented within some inner radial limit 
(Figure 7). In principle, one can use the observed surface-density prohle to estimate the 
amount of mass lost from by the disk throughout its history. However, because large grains 
remain in the disk longer than smaller grains, this estimation implies an implausible amount 
of solid material in the disk. Consequently, either the observed rate of dust generation is 
limited to times near the encounter of the disk with HD 141569A’s binary companions, or 
small grains are created highly preferentially over larger ones. 

Inspired by disk’s large-scale spiral arms, we explored the effect on the disk from a 
parabolic hy-by of the binary companions. Such an encounter is consistent with the mea¬ 
sured proper motion of the system. Our models consider the dynamical evolution of three 
different populations: planetesimals (subject only to gravity), gas (subject to gravity and 
hydrodynamic forces), and dust (subject to gravity, radiation pressure and gas drag). The 
Poynting-Robertson effect, self-gravity, and the dynamical effects of collisions are ignored. 
We vary (3, the initial gas mass, the gas disk opening angle, and the prescription for dust 
generation. We hnd that the dust distribution resembles that of the contemporary plan- 
etesimal distribution. The models are most sensitive to the gas temperature, with hotter 
gas producing more disorganized structures than colder gas. Five-fold differences in the gas 
mass and two-fold differences in (3 do not produce signihcant differences in the simulated 
disks. Quadratic dust generation produces more disorganized and less distinct structures 
than linear dust generation. Smaller initial disks produce smaller and less pronounced spiral 
arms. 

Our favored model (Figure 17) has an initial disk size of 400 AU, /? = 4, and gas with 
with mass 50 and an opening angle of i7/r = 0.1. A planet is introduced to reduce 
the surface density of the planetesimals close to the central star, but the characteristics of 
this planet are not tightly constrained. The model successfully reproduces the large-scale 
structures seen in the scattered-light images of the disk. However, the model does not 
reproduce the tightly-wound spiral rings or other observed asymmetries, which suggests that 
another population of grains (perhaps with smaller j3) or a non-coplanar encounter geometry 
should be investigated. 
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Fig. 1.— Optical depth profile (proportional to the surface density) of the disk shown 
in Figure 5 of 003, derived from the median values of concentric annuli centered on the 
estimated geometric center of the spiral rings. The error bars indicate the standard deviation 
and the dashed lines indicate the upper and lower values in each annulus. 
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Fig, 2,— Scattering opacity at 0,5 /xm as a function of grain size. The scattering opacity is 
proportional to the product the mass opacity (the scattering cross section, o-gca = QscaT^s"^, 
divided by the particle mass, oc s^) weighted by the number of particles (oc and the 

particle mass (oc s^), where s is the particle size. The scattering efficiency Qsca is that of a 
compact astronomical silicate (Laor & Draine 1993; Draine & Lee 1984), To conserve the 
area under the curve with a logarithmic abscissa, the ordinate has been multiplied by s. 
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Heliocentric Velocity (km/sec) 

Fig. 3.— Spectra of J = 2 —1 transition with resolntions of 1 MHz (thick curve) and 

100 KHz (thin curve, after three-point smoothing) from Znckerman et al. (1995). Althongh 
the depth of the central depression is ~ la, it is present in both spectra. 
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Fig. 4.— Radial velocity (divided by the local Keplerian velocity) as a function of radii, for 
/5 = 4 dust particles launched from two different points: Tq = 0.1 (thin lines) and Tq = 0.5 
(thick lines). The behavior of the dust particles with different amounts of gas is shown: no 
gas (highest curve), 10 M®, and 50 M® (lowest). The dotted line is the thermal sound speed 
of the gas. 





- 25 - 


i.O 


0.8 


0.6 

f-i 

0.4 


0.2 

0 


Fig. 5.— Position as a function of time for dust particles released at two different positions: 
0.1 (thin lines) and 0.5 units (thick lines), for a unit disk of gas in thermal equilibrium around 
a 2.3 Mq star. The unit of time is such that the orbital period at the outer edge of the disk 
is 271. (If the disk is 1000 AU in radius, one time unit is 3325 yrs.) The solid lines trace 
particles with /5 = 4, while the dashed lines trace particles with (3 = 10. For each set, three 
lines are shown. From the left: no gas, 10 M 0 and 50 Mq. Notice that for large amounts of 
gas there is little difference between particles with = 10 (the dashed lines) and (3 = A (the 
solid lines). 
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Fig. 6.— Steady-state dust surface density as a function of distance for an isolated disk, 
with dust particles (/5 = 4, s ~ 1 pm according to TAOl or s ~ 3 pm for LL03) for the 
linear model (in which the rate of generation is proportional to the planetesimal snrface 
density). Two sets of snrface densities are shown: one for a disk in which the innermost 
dust generation radius is 0.1 and another for which it is 0.01. The dotted lines indicate the 
surface density with no gas; the dashed lines show the snrface density with 10 Mq of gas; 
the solid lines show the surface density with 50 M 0 of gas. The gas is in LTE. 











r (units of disk radius) 

Fig. 7.— Steady state dust surface density as a function of distance for an isolated disk but 
with the number of dust particles generated proportionally to the square of the planetesimal 
surface density. Line traces as in Figure 6. Notice that a broad ring is produced in all cases. 







Fig. 8.— The time in years it takes for a dust particle inserted with keplerian velocity 
at different radii to reach 500 AU. The central star is 2.3 Mq. The gas sound speed is 
parametrized by the opening angle, H/r, where H is the disk thickness. Solid: /? = 4, gas 
in LTE; dashed: /5 = 4, H/r=0.05; dotted: f3 = 10, gas in LTE. For each set of plots, from 
bottom to top: no gas, 10 M 0 , 100 Mq (this is the gas amount within 500 AU). 
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Fig. 9.— Planetesimal evolution as a function of time, for a prograde encounter. The 
periastron occurs a one unit distance and in this particular example that is also the radius 
of the disk. The units of time are such that the period of a planetesimal at radius one is 
approximately 27r. The closest approach occurs at T=0.0, with the largest of the companions 
closest to the center. The initial surface density distribution goes as The position of 

the binary is indicated with two hlled circles or a dashed line (when outside the presented 
frame). 















Fig. 10.— Gas evolution as a function of time. The gas has H/r=0.1. The parameters of 
the simulation are the same as that of Figure 9. 
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Fig. 11.— The observed profile at a given time is built by adding contributions from dust 
at previous times. In this example, /5 = 4, Mg=5OM0, with H/r=0.1, T = 1.50, quadratic 
dust generation. The panels show the conhguration of dust created at previous times and 
observed at T = 1.50. The plot shows the azimuthally-averaged densities. These plots show 
that the disk is mainly cleared from the inside out. 
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Fig. 12.— Time evolution of the dusty disk. The image shows the predicted density prohle, 
azimuthally smoothed and gaussian convolved to match the measurements. The top row 
shows models with linear dust generation and the bottom row shows models with quadratic 
dust generation. The prohles have been normalized to the same arbitrary maximum value. 
The panels display the logarithm of the density prohles. To highlight faint structure, white 
color in the bottom panels corresponds to pixels 16 times fainter than in the top panels. 
The companions are on top, just outside the frame. The plot inset compares the measured 
density prohle (thin line) with the model density prohle (thick line). To set the distance 
scale, the periastron is set at 718 AU. In these simulations, that is also the size of the initial 
disk. All models have /? = 4, 50 M 0 of gas, H/r = 0.1. Three times are shown: T=0.8, 
T=1.5, and T=2.0. 
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Fig. 13.— The effect of the amount of gas on the dust distribution. Prohles have been 
processed as in Figure 12. Results of different model runs, for T=1.5. All models have 
= 4. The top row shows models with linear dust generation and the bottom row shows 
models with quadratic dust generation. White color in the bottom panels corresponds to 
pixels 4 times fainter than in the top panels. Top row, from left: (1) No gas, linear dust 
generation; (2) 10 M 0 of H/r=0.1 gas, linear dust generation; (3) 50 M 0 of H/r=0.1 gas, 
linear dust generation. Bottom row, from left: (4) No gas, quadratic dust generation; (5) 10 
M 0 of H/r=0.1 gas, quadratic dust generation; (6) 50 M 0 of H/r=0.1 gas, quadratic dust 
generation. In order to correct for artifacts of the dust generation process, the gasless models 
(panels 1 and 4) have been further processed by taking 3-point medians in every pixel. 







Fig. 14.— The effect of the amount of gas on the dust distribution. Same as Figure 13 
but with = 10. The concentric rings for the gasless models are artifacts from the dust 
generation process. 
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Fig. 15.— The effect of the gas temperature on the dust distribution. All models have (3 = A 
and 50 M 0 of gas. The top row shows the results of the linear dust generation method, and 
the bottom row the results of the quadratic one. Top row, from left: (1) H/r=0.05, linear 
dust generation; (2) H/r=0.1 gas, linear dust generation; (3) LTE, linear dust generation. 
Bottom row, from left: (4) H/r=0.05, quadratic dust generation; (5) H/r=0.1 gas, quadratic 
dust generation; (6) LTE, quadratic dust generation. 
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Fig. 16.— Planets and small disks. The top row shows linear dust generation and the bottom 
row shows quadratic dust generation. All simulations have /? = 4, 50 M 0 of H/r = 0.1 gas. 
(1) 400 AU initial disk, linear dust generation; (2) 400 AU initial disk, linear dust generation, 
with an eccentric (e = 0.6), 5 Mj^p planet with semimajor axis 100 AU; (3) 400 AU initial 
disk, quadratic dust generation; (4) 400 AU initial disk, quadratic dust generation, with an 
eccentric 5 Vijup planet within 160 AU; 
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Fig. 17.— Comparison between best model and observations. Both are logarithmic stretches. 
The model {{3 = 4, 50 M 0 of H/r = 0.1 gas, cuadratic dust generation, with an eccentric 5 
Mjup planet within 100 AU) has been rotated counterclockwise by 13 degrees. 
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Fig. 18.— Comparison between measured and modeled profiles, at T = 1.50. The traces 
show density prohles generated quadratically in small disks {j3 = 4, 50 M 0 of H/r = 0.1 
gas), with and without planets (Figure 16, panels 3 and 4). Also shown are the planetesimal 
and gas surface densities. These are steeper than initially set because of the truncating effect 
of the encounter. 










